Anomalous exponents in the rapid-change model of the 
passive scalar advection in the order e 3 . 
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Field theoretic renormalization group is applied to the Kraichnan model of a passive scalar 
advected by the Gaussian velocity field with the covariance (v(t, x)v(t ,x)) — (v(t, x)v(i', x')) oc 



Inertial-range anomalous exponents, related to the scaling dimensions of tensor 



composite operators built of the scalar gradients, are calculated to the order e of the e expansion. 
The nature and the convergence of the e expansion in the models of turbulence are briefly discussed. 
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The investigation of intermittency and anomalous scal- 
ing in fully developed turbulence remains essentially an 
open theoretical problem. Both the natural and numeri- 
cal experiments suggest that the deviation from the pre- 
dictions of the classical Kolmogorov theory Q is even 
more strongly pronounced for a passively advected scalar 
field than for the velocity field itself; see, e.g., Ref. ||] 
and literature cited therein. At the same time, the prob- 
lem of passive advection appears to be easier tractable 
theoretically: even simplified models describing the ad- 
vection by a "synthetic" velocity field with a given Gaus- 
sian statistics reproduce many of the anomalous features 
of genuine turbulent heat or mass transport observed in 
experiments. Therefore, the problem of passive scalar 
advection, being of practical importance in itself, may 
also be viewed as a starting point in studying anomalous 
scaling in the turbulence on the whole. 

Most progress has been achieved for the so-called 
rapid-change model ||: for the first time, the anoma- 
lous exponents have been calculated on the basis of a 
microscopic model and within regular perturbation ex- 
pansions; see Refs. p|-|r5| and references therein. 

In that model, the advection of a passive scalar field 
9(x) = 0(t, x) is described by the stochastic equation 



d t 9 + (Vidi)6 = voAO + /, 



(1) 



where dt = d/dt, di = d/dxi, vq is the molecular diffusiv- 
ity coefficient, A is the Laplace operator, v(x) = {vi(x)} 
is the transverse (owing to the incompressibility) velocity 
field, and / = f(x) is an artificial Gaussian scalar noise 
with zero mean and correlator 



(f(x)f(x')) = S(t-t')C(r/L), 



r = x — x 



(2) 



The parameter L is an integral scale related to the scalar 
noise, and C(r/L) is some function finite as L — > oo. 
Without loss of generality, we take L = oo and C(0) = 1. 

In the real problem, the field v(x) satisfies the Navier- 
Stokes equation. In the rapid-change model it obeys a 
Gaussian distribution with zero mean and correlator 



{v l {x)v j (x')) = D 



6(t - 1') 



(2 

xk- d 



exp[ik • (x — x )], 



(3) 



where -fV,(k) 



kikj /k 2 is the transverse projector, 



k = |k|, Do > is an amplitude factor, d is the dimen- 
sionality of the x space, and < e < 2 is a parameter 
with the real ("Kolmogorov") value e — 4/3. The in- 
frared (IR) regularization is provided by the cut-off in 
the integral ([|) from below at k ~ m, where m = \ ji is 
the reciprocal of another integral scale £; the precise form 
of the cut-off is not essential. The relation Dq/vq = A £ 
defines the characteristic ultraviolet momentum scale A. 

The issue of interest is, in particular, the behavior of 
the equal-time structure functions 



S n (r) = (l9(t,x)-8(t,x')r), r=|> 



(4) 



in the inertial-convective range A»l/r»m. 

In the isotropic model (0)-(||), the odd multipoint cor- 
relation functions of the scalar field vanish, while the 
even equal-time functions satisfy linear partial differen- 
tial equations [||-|5| . The solution for the pair correlator is 
obtained explicitly; it shows that the structure function 
S2 is finite for m — jjj. The higher-order correlators 
are not found explicitly, but their asymptotic behavior for 
to — * can be extracted from the analysis of the nontriv- 
ial zero modes of the corresponding differential operators 
in the limits l/d -> f§, e -> [|§, or s -> 2 §,§. 
It was shown that the structure functions in the inertial- 
convective range exhibit anomalous scaling behavior: 



S 2 »(r)ocIV^" (2 ~ e) (mr) 1 



(5) 



with negative anomalous exponents A„ , whose first terms 
of the expansion in l/d and e || have the forms 



-n(n~2)e/2d+0(l/d 2 ) = 
-n{n ~ 2)e/2{d + 2) + (9(e 2 ). 



(6) 



In paper the field theoretic renormalization group 
(RG) and operator product expansion (OPE) were ap- 
plied to the model dl|)-(p|). In the RG approach, the 
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anomalous scaling for the structure functions and vari- 
ous pair correlators is established as a consequence of the 
existence in the corresponding operator product expan- 
sions of "dangerous" composite operators (powers of the 
local dissipation rate) , whose negative critical dimensions 
determine the anomalous exponents A„ . The anomalous 
exponents were calculated in Ref. 0] in the order e 2 of 
the e expansion for the arbitrary value of d; generaliza- 
tion to the compressible case was given in 10 111]. The 
main advantage of the RG approach (apart from its cal- 
culational efficiency) is the universality: it can equally be 
applied to the case of finite correlation time . 

In this paper, we present the anomalous exponents 
and other quantities for the model (0)-(§) in the order 
e 3 . Here we give only basic ideas and results; more ex- 
haustive discussion of the calculational technique will be 
given elsewhere. A general review of the RG approach 
to the statistical theory of turbulence can be found in 
Rcfs. [ ^6|Jl7| l; the case of the Kraichnan model is discussed 
in (|] in detail. 

The stochastic problem ([[])- (||) can be reformulated as 
a multiplicatively renormalizable field theoretic model; 
the corresponding RG equations have an IR attractive 
fixed point. This implies existence of the infrared scaling 
behavior for all correlation functions with certain scaling 
dimensions, calculated as series in e (in this sense, the 
exponent e plays in the RG approach the same role as 
the parameter e = 4 — d does in the RG theory of critical 
behavior). In particular, for the structure functions (Q), 
(|J) in the IR asymptotic range (Ar 3> 1) one obtains 

^ (2 ~ e) Xn(mr). 



S 2n (r) cx iV r " 



(7) 



The behavior of the scaling functions Xn (mr) at mr 
(inertial-convective range) is obtained with the aid of the 
operator product expansion: 



Xn (mr) 



F 



(mr) 1 



(8) 



where the sum runs over all possible composite opera- 
tors F entering the OPE for a given structure function, 
Ap are their critical dimensions, and Cf are numerical 
coefficients analytical in (mr) 2 and finite at mr = 0. 

The key role is played by the critical dimensions A„;, 
associated with the tensor composite operators 



d tl e ■ ■ ■ d^e (dmt 



(9) 



where / is the number of the free vector indices and 
n = I + 2p is the total number of the fields 6 entering 
the operator; the vector indices of the symbol F n i are 
omitted. 

The dimension A„ = A„o of the scalar operator is 
nothing other than the anomalous exponent in Eq. (|^); 
see Ref. ||. The dimensions with / ^ become relevant 
if the forcing ^ is anisotropic: A n i corresponds to the 
zero-mode contribution to the l-th term of the Legendre 
decomposition for the function S n ; see Ref. They 
can be systematically calculated as series in e: 



(10) 



k=l 



with the first-order coefficient 111 31 



nl 



(n-2) (d+l)l(d + l -2) 



2(d + 2) 2(d-l)(d+2) 



(11) 



(see also Refs. [y_3|jlj]; for I = this gives the result of 
||, while for n — 3 and I — 1 or 3 the result of is 
recovered) . 



The coefficients A^ 2 J and were obtained in Ref. |9[] 
for any n and d; the result for general I is presented in 
|ll| . In particular, one has 



(2) 



(2) 



n(n - 2)(0.000203n - 0.02976) 
J 2 (0.01732n + 0.01223) 



(12) 



for d = 2 and 

a ( 2 ) = n (n- 2)(0.00203n- 0.00384) 
- 1(1 + l)(0.00710n - 0.00619) 



nl 



(13) 



for d = 3 (analytical results are too cumbersome and will 
not be given here; see Refs. pjLl||). 

Now let us turn to the 0(e 3 ) contribution. No analyt- 
ical result for it is available for general d; the numerical 
results have the forms 



n(n - 2)(0.005472n 2 + 0.0649n + 0.0647) + ; 2 (-0.02161n 2 - 0.1023n + 0.2406 + 0.01841/ 2 ) 



for d = 2 and 

A ni = n ( n - 2)(0.00140ti 2 + 0.0199n + 0.0343) + l(l + l)(- 0.00420n 2 - 0.0241n + 0.0028(Z 2 + I + 18)) 



(14) 



(15) 



for d = 3. The quantity A^ ; can be expanded as a series 
in 1/d; the coefficients of such expansion can be found, 



in principle, to any given order. For general n and I to 
the order 1/d 2 we have obtained 

A„ ; = e[ - n(n - 2)(1 - 2/d)/2d+ (Z/2)(l - 2/d + l/d + 2/d 2 )] + 3e 2 (n - 2)(n - I) /Ad 2 + 
+ e 3 (n - I) [l.74988(n - 2) - 0.624916Z]/d 2 . 



(16) 
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Note that the £ 2 and £ 3 contributions decay for d — ► oo 
faster than l/d in agreement with the 0(1/ d) result ob- 
tained in Ref. B for A„q. Moreover, from Eq. (16) it 



follows that the leading 0(1/ d 2 ) terms in these contri- 
butions vanish for n = I, so that the decay at d — > oo 
becomes even faster: 



en 



/2 + n(n - 1) {e /(d - l)(d + 2) - e 2 [l + (2n - 7)/d] /d 3 - £ 3 (3n - 8)/2d 4 } + 0(e 4 ), 



(17) 



with the accuracy of 0(l/d 4 ). 

We also recall that A20 = to all orders in e in agree- 
ment with the exact solution for the second-order struc- 
ture function || , and that the exact nonperturbative re- 
sult for A22 exists for all e and d Jlj. 

For the isotropic model only scalar operators 

enter the expansion (||) , the number of the fields 9 in the 
operators does not exceed the number of 0's on the left 
hand side, and the leading term of the small-mr behavior 
is given by the operator with the minimal dimension Ap- 
This allows one to identify the anomalous exponent A n 
in Eqs. (||) and (^|) with the critical dimension A„o of 
the scalar operator F n Q. 

If the noise covariance (|J) involves some fixed constant 
vector n (large-scale anisotropy), the above results for the 
dimensions A n i do not change, but the operators with 
I =^ also enter the right hand side of Eq. (||) and give 
rise to contributions proportional to Pi(z), the Z-th order 
Legendre polynomial, z being the angle between the vec- 
tors n and r. The odd structure functions S^n+i become 
nontrivial, and the leading term of their inertial-range 
behavior is determined by the dimension A2 n +i,i of the 
vector operator F 2n +iA 

S 2n +i{r) oc D- ,l - 1/2 r («+V2)(2-e) ( mr )A 2 „ +1 ,i (18) 

(for more detail, see Refs. fl2]| ). 

In Fig. 1, we show the dimension A4 (which determines 
the anomalous exponent for S4) for d = 3 in the first, 
second, and third orders in e. In Figs. 2 and 3, we show 
the "anomaly" 7, defined by the relation S3 oc r 3 r , for 
d = 3 and 2, respectively; note that the 0(e 2 ) curve lies 
above the Oie) line for n — 3 and below it for n = 4. In 
the same Figures, we also present nonperturbative results 
obtained for n — 4 in Refs. |^| using numerical simula- 
tions, and for n = 3 in |Q using numerical integration of 
the zero-mode equations (A4 = (4 — 2(2 in the notation 
of M and 7 = 3 — A in the notation of [Q). 

An important issue which can be discussed on the ex- 
ample of the rapid-change model is that of the nature and 
convergence properties of the e expansions in models of 
turbulence and the possibility of their extrapolation to 
finite values e ~ 1. Figures 1-3 show that the agreement 
between the e expansion and nonperturbative results for 
small e improves when the higher-order terms are taken 
into account, but the deviation becomes remarkable for 
e ~ 1 and decreasing d. Furthermore, the coefficients of 
the e series appear more irregular for d = 2 (see Fig. 3), 
while the forms of the nonperturbative results are 



not much affected by the choice of d. 

Such behavior can be understood on the basis of the 
exact analytical result for A22, which can be written in 
the form Q 

2A 22 = -d - 2 + £ + y/(e + £+)(£ + £_), (19) 
where 

e± = (d 2 + d + 2 ± y/8d(d+l))/(d - 1). 

It shows that the corresponding e expansion has the fi- 
nite radius of convergence £_, ranging from to 00 when 
d varies from 1 to 00; in particular, £_ ~ 1.1 for d = 2 
and £_ ~ 2.1 for d = 3. Hence the naive summation of 
the £ expansion for A22 works only in the interval e < £_ , 
which decreases almost linearly with (d — 1). In order to 
recover the behavior of A22 from its e series for larger e, 
it is necessary to isolate explicitly the singularity at £_ in 
Eq. (|l9]), thus changing to a kind of improved £ expan- 
sion (whose radius of convergence becomes £+ 3> £-). 
In practice, the first three terms of this improved expan- 
sion approximate the exact result ( |l9|) equally good for 
all < £ < 2, both in two and three dimensions. 

The difference with the models of critical phenomena, 
where e series are always asymptotical, can be traced 
back to the fact that in the rapid-change models, there is 
no factorial growth of the number of diagrams in higher 
orders of the perturbation theory. The divergence for 
c? — s> 1 is naturally explained by the fact that the trans- 
verse vector field does not exist in one dimension (we 
also recall that the RG fixed point diverges at d = 1; 
see Ref. Thus it is natural to assume that the se- 

ries for higher-order exponents A n i also have finite radii 
of convergence with the behavior similar to that of £_. 
Therefore, in order to obtain reasonable predictions for 
finite values of £, one should augment plain e expansions 
by the information about the location and character of 
the singularities. Such information can be extracted from 

(k) 

the asymptotical behavior of the coefficients A^ ; in Eq. 
© at large k. To our knowledge, this problem has never 
been studied for dynamical models like (|l])-(|3|); the in- 
stanton analysis developed in Refs. jl5| has mostly been 
concentrated on the behavior of the exponents in the 
limit n — > 00. One can hope that the implementation 
of the instanton calculus within the RG framework will 
give the solution of this important problem. 
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FIG. 1. The dimension A4 for d — 3 vs e: the O(e), 
0(e 2 ) and 0(e 3 ) approximations (from above to below). 
Dashed line: numerical simulation by Refs. [8]. 
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FIG. 2. The exponent 7 for d = 3 vs e: the 0(e 2 ), 
O(e) and 0(e 3 ) approximations (from above to below). 
Dashed line: numerical solution by Refs. [7]. 
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FIG. 3. The exponent 7 for d — 2 vs e: the 0(e 2 ), 
0(e) and 0(e 3 ) approximations (from above to below). 
Dashed line: numerical solution by Refs. [7]. 
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